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Bose-Einstein condensates of 10 4 85 Rb atoms in a cylindrical trap are studied using a recently 
proposed modified Gross-Pitaevskii equation. The existence of a Feshbach resonance allows for 
widely tuning the scattering length of the atoms, and values of the peak gas parameter, x p k, of 
the order of 10 -2 can be attained. We find large differences between the results of the modified 
Gross-Pitaevskii and of the standard Thomas-Fermi, and Gross-Pitaevskii equations in this region. 
The column densities at z — may differ by as much as ~ 30% and the half maximum radius 
by ~ 20%. The scattering lengths estimated by fitting the half maximum radius within different 
approaches can differ by ~ 40%. 
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Bose-Einstein condensation (BEC) of magnetically trapped alkali atoms has been achieved in several experimental 
setups, most of them in regimes where the atomic gas is considered to be very dilute, i.e. the average interatomic 
distance is much larger than the range of the interaction. As a consequence, the physics is dominated by two body 
collisions, generally well described in terms of the s-wave scattering length, a. The crucial parameter defining the 
condition of diluteness is the gas parameter, x(r) — n(r)a 3 , where n(r) is the local density. For low values of the 
average gas parameter, x av < 10~ 3 , the mean field Gross-Pitaevskii (GP) equation |Q is the logical tool to study the 
system. 

The gas parameter can be brought outside the regime of validity of the GP equation by two ways: either by 
increasing the number of atoms, N, in the condensate, or by changing their effective size. Recent experiments have 
explored both possibilities. On one side they have reached very large N values, N ~ 10 8 ; on the other the scattering 
lengths have been widely tuned. The second approach promises to be much more efficient to reach large x regions. 

In a recent experiment performed at JILA (^H, it was possible to confine about 10 4 atoms of 85 Rb in a cylindrical 
trap. By exploiting the presence of a Feshbach resonance at a magnetic field of B ~ 155 Gauss, the scattering length 
was varied from negative to very high positive values. Actually, it was suggested some time ago || that a could 
be modulated by taking advantage of its expected strong variation in the vicinity of a magnetic induced Feshbach 
resonance in collisions between cold alkali atoms. Several experiments have supported this proposal by demostrating 
this type of scattering length variation for several alkali atoms, as 85 Rb, Cs and Na ■ However, only very recently 

the Boulder group at JILA has been succesful in producing stable 85 Rb condensates where a can be effectively tuned 
over a very wide range. 

To explore with confidence such regimes of high values of the parameter x, a necessary task is to investigate the 
accuracy of the GP equation. Moreover, to quantify the limitations of the GP description is of particular relevance 
since the empirical estimate of the scattering length is based upon the so called Thomas-Fermi (TF) approximation 
to the GP equation. The TF approximation amounts to disregard the kinetic energy term in the GP equation. The 
TF and GP results are expected to coincide for large N and/or a values. In the experimental analysis, a common 
procedure consists in measuring the column density, given by the integral of the particle density along a direction 
perpendicular to the simmetry axis of the trap, n c (z) — J dx n(x, 0, z). The x direction coincides with that of the 
light beam used to image the atomic cloud. Then, the scattering length is inferred by finding the value of a that, in 
the framework of the TF equation, provides a column density with the same experimental size. 

In this paper, we will use a recently proposed modified Gross Pitaevskii equation (MGP) (8) to estimate the 
corrections to the GP results. We begin by briefly recalling the derivation of the MGP equation. Then, in the simpler 
case of a spherical trap we will compare different approaches in a situation where the number of atoms and the 
frequency of the trap are kept fixed and the scattering length is allowed to vary in a representative range of values. 
This study is done to ascertain the degree of reliability of the MGP results. Finally, we will consider a cylindrical 
trap corresponding to the experimental situation. Depending on the value of the scattering length, the corrections to 
the GP results can be as large as 30 % for energy, chemical potential and ~40 % for the extracted scattering length. 

In the spherical case, the energy functional associated with the MGP equation: 
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is obtained in the local-density approximation by keeping the first two terms in the low density expansion for the 
energy density of a homogeneous system of hard-spheres, whose diameter coincides with the scattering length 0: 
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Up to this order of the expansion, the details of the potential do not show up, and any potential with the same 
scattering length would give identical results. This universal behavior has recently been checked by a diffusion Monte 
Carlo calculation (DMC) ]l0| ], providing the exact solution of the many-body Schrodinger equation. In Refs. |p|,pO 11 
it was shown that, for an uniform system, the first term of the expansion is accurate only at very low values of x, 
while the addition of the second term gives a good representation of the exact DMC results up to x — 1CP 2 . The 
inclusion of the logarithmic term severely spoils the agreement already at intermediate x-values and therefore it has 
not been incorporated into the functional energy, Emgp W\ ■ ^ ^ ne same references it was also shown that the energy 
functional computed in hypernetted chain (HNC) theory, explicitly taking into account the interatomic correlations 
induced by the potential, provided a description very close to the MGP one. 

It is convenient to simplify the notation by expressing lengths and energies in harmonic oscillator (HO) units. 
The spatial coordinates, the energy, and the wave functions are rescaled as r = a£ro r i) E — huEi, and \&(r) = 
(iV/a| ?0 ) 1//2 \E'i(ri), where 'l'i(ri) is normalized to unity and ano = (h/mto) 1 ^ 2 . Using these new variables and 
performing a functional variation of EmgpW\i one S e ^ s the modified Gross-Pitaevskii equation, 
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where a\ — a/auo and £ii is the chemical potential in HO units. The CP approximation is recovered by dropping 
the | ipi(ri) | 3 term. 

Table | gives some results for N — 10 4 85 Rb atoms confined in a spherical trap with an oscillator angular frequency 
uho/Zk = ( w l w z) 1,/3 /27r = 12.77 Hz, where uj±/2n = 17.5Hz and u> z /2tt = 6.9 Hz are the radial and axial frequencies 
associated with the external potential of the cylindrical trap used in Ref. (|] . In the table we study the dependence of 
the energy per particle and of the chemical potential on the scattering length, given in units of the Bohr radius of the 
Hydrogen atom, ay. For this trap, ciho — 57657ao. We also show the TF results. In this approximation it is often 
possible to derive simple analytical expressions p"2| , useful to get quick estimates of several quantities. For instance, 
[iJ F = l/2(15aiiV) 2 / 5 . Also reported are the peak values of the gas parameter, x p t — n.(0) a 3 = Na\ | V>i(0) | 2 , whose 
TF estimate is x^ = (15 2 a} 2 A^ 2 ) 1 / 5 /(87r). 

At low values of the scattering length, MGP corrections are small and the TF approach to the GP equation is not 
fully satisfactory. As expected, the TF and GP results are much closer when a increases and the MGP corrections 
become important and of the order of 30% at the largest value of a — lOOOOao. x p k increases with a. The MGP density 
distribution gets wider and x^| GP is depleted with respect to both TF and GP because of the repulsive character of 
the extra term in the MGP equation. The HNC results are confortably close to the MGP ones, supporting the use of 
only the latter approach in the remaininng of the paper. 

An analogous behavior is found in the anisotropic case. The results for the cylindric trap used in the experiments 
are given in Table H Most of the TF results are again analytical. The HO units in the cylindric case are: r = 
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(H/muj±, ) 1 / 2 . The trap deformation parameter is 



a± tH c>ri, E = Twj±Ei, *(r) 
A = u z /u± = 0.39. 

An accessible experimental quantity connected to the density profile is the already defined column density. Its TF 
expression is 
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A measure of the extension of the condensate is the half maximum radius of the column density, R\ /2 , defined as the 
z\ value where n c (z\ = R1/2) = \n c (Q). Also interesting is the full strength at half maximum, FSHM, given by the 
integrated strength of the column between the ±i£i/2 values. 

In Fig. ^, we show the column densities in different approaches, for the same set of scattering lengths reported in the 
Tables. The solid and dashed lines correspond to the MGP and TF results, respectively; stars are the GP densities; 
the triangles in the two upper panels give n c (z%) evaluated in MGP, but changing the scattering length to reproduce 
R1/2.TF1 supposedly corresponding to the measured radius. The two values are a/ao — 5920 for R 1 / 2 ,tf=10-20 and 
a/ao = 4940 for i?i/2,i\F=9-75 and the related MGP columns are practically identical to the TF ones. 
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The GP and TF results almost coincide and the MGP corrections are more sizeable at the two largest values of a, 
where x p k becomes of the order of 10~ 2 . Because of the repulsive nature of the MGP extra term, R\/2,mgp is larger 
than R\/2,tf, while FSHMmgp is smaller than FSHMtf, and for low zi-values n^ IGP (zi) lies below n^ F (zi). A 
smaller scattering length is required to reproduce R\/2,tf an d FSHMyp in the MGP approach in the high x p k region. 
In fact, we find a reduction of ~ 40% of a/a,Q. This analysis shows that using the TF column density to extract the 
scattering length in the large gas parameter regime could lead to severe overestimates in this kind of trap geometry. 

The lower panels of the figure roughly correspond to x p k ~H0~ 3 -10~ 4 . As expected, the MGP corrections are 
smaller and the computed R1/2 values become closer when a/ao decreases. 

Fig. |^ shows the scattering length as a function of FSHM and R1/2 for the cylindrical trap within the three 
methods we have analyzed. The figure stresses that, depending on the FSHM and R1/2 values and on the approach, 
the estimates of a can differ by up to 40%. 

In conclusion, we find that the MGP equation induces corrections of 30% in the ground state properties of the 
condensate, when the conditions of the JILA experiments for 85 Rb are considered. Comparable corrections are 
obtained for the column densities, where large differences between the MGP and the standard TF and GP results 
may be found. These differences appear to be relevant for the extraction of the scattering length when large values 
of the gas parameter come into play. MGP is still a mean field theory, since it tries to incorporate correlation effects 
into the average single particle potential. However, we believe that its predictions are probably indicative at those 
regimes attained in recent experiments. 
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TABLE I. Ground state properties of N = 10 4 85 Rb atoms confined in a spherical trap (o>/27r=12.77 Hz) in different 
approaches. /ii=chemical potential, Ei/N= energy per atom, :r p fc=peak gas parameter. Energies in HO units. 



a/ao 


1400 


3000 


8000 


10000 


mP' 
M? P 

,.MGP 
Ml 

HNC 
Ml 


13.29 
13.41 
13.95 
13.90 


18.02 
18.12 
19.82 
19.66 


26.00 
26.75 
33.34 
33.33 


29.18 
29.24 
38.01 
38.41 


Ei /N 
Ef GP /N 
E? NC /N 


9.50 
9.66 
10.00 
9.97 


12.87 
13.01 
14.09 
13.98 


18.57 
19.16 
23.40 
23.24 


20.84 
20.93 
26.60 
26.59 


^pk 

%pk 
„MGP 
•^pk 

x pk 


6.23 xl0~ 4 
6.26 xl0~ 4 
5.70 xl0~ 4 
5.75 xl0~ 4 


3.88 xl0~ 3 

3.89 xl0~ 3 
3.18 xl0~ 3 
3.24 xl0~ 3 


4.09 xl0~ 2 
4.09 xl0~ 2 
2.59 xl0~ 2 
2.52 xlO" 2 


6.98 xlO -2 

6.99 xl0~ 2 
4.10 xl0~ 2 
3.85 xl0~ 2 



TABLE II. Ground state properties of N = 10 4 85 Rb atoms in the cylindrical trap described in the paper. Energies in HO 
units. 



a/a 


1400 


3000 


8000 


10000 


M? P 

MGP 

Mi 


9.70 
9.82 
10.22 


13.15 
13.25 
14.51 


19.47 
19.55 
24.38 


21.29 
21.36 
27.79 


Ei"'/N 
E GP /N 
Ef IGP /N 


6.93 
7.08 
7.33 


9.39 
9.52 
10.31 


13.91 
14.00 
17.09 


15.21 
15.29 
19.43 


^pk 
MGP 
x pk 


6.23 xl0~ 4 
6.28 xl0~ 4 
5.72 xl0~ 4 


3.88 xl0~ 3 
3.90 xl0~ 3 
3.19 xl0~ 3 


4.09 xl0~ 2 

4.10 xl0~ 2 
2.60 xl0~ 2 


6.98 xlO" 2 
7.00 xl0~ 2 
4.10 xl0~ 2 
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FIG. 1. Column densities at four values of the scattering length for the cylindrical trap. Dashed lines= TF, stars= GP, solid 
lines= MGP. The triangles in the first (second) upper panel give the MGP column density at a/ao=5920 (4940). 
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